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Statistical mechanics and large-scale velocity fluctuations of turbulence 

Hideaki Mouri, a ) Akihiro Hori, b ) Yoshihide Kawashima, b ) and Kosuke Hashimoto* 5 ) 

Meteorological Research Institute, Nagamine, Tsukuba 305-0052, Japan 

(Dated: 18 November 2011) 

Turbulence exhibits significant velocity fluctuations even if the scale is much larger than the scale of the 
energy supply. Since any spatial correlation is negligible, these large-scale fluctuations have many degrees of 
freedom and are thereby analogous to thermal fluctuations studied in the statistical mechanics. By using this 
analogy, we describe the large-scale fluctuations of turbulence in a formalism that has the same mathematical 
structure as used for canonical ensembles in the statistical mechanics. The formalism yields a universal law 
for the energy distribution of the fluctuations, which is confirmed with experiments of a variety of turbulent 
flows. Thus, through the large-scale fluctuations, turbulence is related to the statistical mechanics. 



I. INTRODUCTION 

Turbulence is induced by supplying kinetic energy at 
some scale L. This energy could be transferred to both 
the larger and the smaller scales. 1 ' 2 However, as sketched 
in Fig. 1, the energy is on average transferred to smaller 
and smaller scales because it is eventually dissipated into 
heat at the smallest scale, i.e., the Kolmogorov length r\. 
The energy transfer from L to rj consists of many random 
steps, each of which occurs preferentially between neigh- 
boring scales. 1-3 Hence, although motions at the scale L 
depend on the flow configuration for the energy supply, 
such dependence is lost during the energy transfer. The 
resultant small-scale motions exhibit universal features, 
which have been studied as a representative of spatially 
correlated nonequilibrium fluctuations. 3 

The kinetic energy could be transferred to scales much 
larger than L and could cause velocity fluctuations there 
(Fig. 1). Also in these large-scale fluctuations, we expect 
some universality. To lose dependence on the flow config- 
uration, the energy transfer could have a sufficient num- 
ber of random steps. Any step prefers to occur between 
neighboring scales because such scales alone interact in a 
coherent manner. 3 

If the turbulence is stationary and is homogeneous at 
least along one direction, we expect that the large-scale 
fluctuations are analogous to thermal fluctuations in an 
equilibrium state described in the standard textbooks 4-6 
of the statistical mechanics. The stationarity implies that 
no net energy is transferred across the large scales, while 
the homogeneity implies that no net energy is transferred 
in space along that direction. This is analogous to the 
case of the thermal equilibrium at constant temperature, 
where occurs no net heat transfer. 4-6 In addition, at the 
large scales of the homogeneous turbulence, we could ig- 
nore any of the spatial correlations. Then, the fluctuation 
energy is additive. Its value for a large-scale region is the 
sum of its values for the yet large parts of the region that 
are not correlated at all. 4 ' 5 The large-scale fluctuations 



are thus considered as a collection of many distinct mo- 
tions. This is again analogous to the case of the thermal 
fluctuations, which have many degrees of freedom. 

The large-scale fluctuations of turbulence are known to 
be significant, regardless of the flow configuration, 7-10 as 
considered by Landau. 3,10-12 However, their details are 
still not known. Experimentally or numerically, any de- 
tailed study needs long data for many realizations of the 
large scales. Such data have not been available. The situ- 
ation is nevertheless improving, 9 owing to improvements 
in experimental technologies. 

By assuming the universality and the additivity, we de- 
scribe the equilibrium large-scale fluctuations of the sta- 
tionary homogeneous turbulence in a thermostatistical 
formalism, i.e., a formalism that has the same mathemat- 
ical structure as used for the statistical mechanics. 13-15 
The formalism is confirmed with long experimental data 
of a variety of turbulent flows obtained in a wind tun- 
nel. We thereby demonstrate that turbulence is related 
through its large-scale fluctuations to the statistical me- 
chanics. 



II. CONFIGURATION AND COARSE GRAINING 

Let us consider a lateral velocity v(x) along some one- 
dimensional cut x of three-dimensional stationary turbu- 
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FIG. 1. Sketch of three-dimensional stationary turbulence. 
The open arrows indicate the local and instantaneous energy 
transfer. The filled arrows indicate the net energy transfer 
along with the energy supply and the energy dissipation. 
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lence. The longitudinal velocity u(x) is to be also consid- 
ered, by replacing v with u in the following descriptions. 
The turbulence is assumed to be homogeneous along the 
one-dimensional cut. The average (v) is subtracted so as 
to have (v) = anywhere below. As a characteristic scale 
L of the energy supply, we use the correlation length of 
the local energy v 2 . The usual definition is 



L„2 — 



_ ^{[v 2 {x + r)-{v 2 )][v 2 {x)-{v 2 )])dr 



but our definition for later convenience is 

_ ((v 2 (v 2 )) 2 ) - 

L - 2 - ^2)2 



(lb) 



We have L v 2 = L v 2 in a special case where the distribu- 
tion of v is Gaussian, (v 4 ) = 3(v 2 ) 2 . 

The one-dimensional cut is divided into segments with 
length R. For each segment, the center of which is ten- 
tatively defined as x c , the energy v 2 is averaged as 



I r+R/2 
v r( x c) = s / v 2 (x c +x)dx. 

R J-R/2 



(2) 



We focus on this coarse-grained energy. The mean square 
of v R around its average (v R ) = (v 2 ) is 



4,16 



<K-<4» 2 > 

R 



(3a) 



= &j a (R-r)([v 2 (x + r)-(v 2 )][v 2 (x)-(v 2 )])dr, 

where the same symbol (•) is used to denote averages over 
the positions x and over the segments. Since the energy 
is supplied at around the scale L v 2, we assume that any 
n-point spatial correlation of v 2 decays fast enough to 
become negligible at r >• L v 2. In other words, we assume 
additivity of Rv 2 R at R > L v 2. Then, Eqs. (1) and (3a) 
yield 



<(4-<4» 2 > 



4L„ 



R 



^-(v 2 ) 2 at R^>L V 2. (3b) 



The assumption also implies ((v R — (v R )) n ) oc (v )" for 
n = 3, 4, where ((v R — (v R )) n ) in itself depends on 
the spatial correlations of v 2 among up to n points. 16 We 
specify the coefficients of these relations with a thermo- 
statistical formalism [Eq. (10)]. 



III. THERMOSTATISTICAL FORMALISM 

There is an analogue of Eq. (3b) in the statistical me- 
chanics of equilibrium systems with many degrees of free- 
dom. It is a formula for thermal fluctuations of the energy 
£ in a canonical ensemble that has the size R and is in 
contact with a heat bath at temperature T: 4,6 



((E - (E)) 2 ) = C R T 2 with d 



The derivative is taken for constant R. We have assumed 
that Rv 2 R at R » L v 2 is additive. Since E is also additive, 
Eq. (4) is equivalent to Eq. (3b) through the correspon- 
dences 



T = 



and E = N 



-(l-^C)(v 2 )], (5a 



with 



N = 



R 

4L„: 



»1, 



and hence 



(E) = (NT and C R = (N. 



(5b) 



(5c) 



Here ( > is an arbitrary constant that is to be deter- 
mined later [Eq. (8)]. 

Each segment with length R consists of N subsegments 
with length AL V 2 and mean energy y/C,{v 2 ). The adjacent 
subsegments might be somewhat correlated, but such a 
correlation has to be negligible at the larger scales. Thus, 
the segment is a collection of N distinct motions, which 
are individually attributable to the energy-containing ed- 
dies. Once determined, N = R/AL V 2 is kept constant 
even if R varies afterwards [Eq. (12d)], by assuming that 
the turbulence expands or contracts in a self-similar man- 
ner so that L v 2 varies with R. The segment is in an equi- 
librium with the surrounding turbulence that serves as a 
heat bath at T = (v 2 )/y/(. Although this is not a true 
temperature, the analogy is close enough to reproduce 
the observed distribution of v 2 R in Sec. IV. 

The energy distribution P(E) in any canonical ensem- 
ble is determined by the heat capacity C R , 4 ' 6 through a 
series of basic relations of the statistical mechanics. Since 
Cr is related to the entropy (S) as Cr = T(dT(S))n, we 
integrate Cr = (N in Eq. (5c) to obtain 



(S) = (N 



In 



1 



(6a) 



with a constant of integration (N(l — In To) that could 
depend on R via To. The Hclmholtz free energy (F) = 
(E)-T(S) is 



(F) = -CNTln 



The partition function Z = exp(— (F)/T) is 



T 



Civ 



(6b) 



(6c) 



From the inverse of the Laplace transformation Z(T) = 
J °° Q(E) exp(—E/T)dE, we obtain the density of states 

Sl(E) = ^ Ar - 1 /r(C^V)T fAr , where T is the Gamma func- 
tion. Lastly, P(E) = Sl(E)exp(-E/T)/Z(T) is obtained 
independently of T as 



P(E) = 



E< N - 1 exp(-E/T) 
T((N)TC N ' 



(7) 
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TABLE I. Experimental conditions and turbulence parameters of grid turbulence (Gl, G2, and G3), boundary layer (Bl, B2, and 
B3), and jet (Jl, J2, and J3). The velocity derivative was obtained as d x v — [8v(x+5x)—8v(x—Sx)—v(x+26x)+v(x—26x)]/125x 
with the sampling interval 8x = U/f. 



Quantity 






Units 


C 1 




bo 


Bl 


B2 


T=> O 
-DO 


J 1 


T9 
J Z 


J o 


Measurement position 






m 


+ 1.5 


+ 1.5 


+2.0 


+ 12.5 


+ 12.5 


+ 12.5 


+ 15.5 


+ 15.5 


+ 15.5 


Measurement position 






m 


1.00 


1.00 


1.00 


0.35 


0.30 


0.25 


0.40 


0.40 


0.40 


Sampling frequency 


f 




kHz 


10 


24 


54 


10 


26 


60 


16 


44 


90 


Total number of data 






10 s 


1 


4 


4 


1 


4 


4 


1 


4 


4 


Kinematic viscosity 


V 




cm 2 s _1 


0.141 


0.142 


0.143 


0.138 


0.143 


0.143 


0.139 


0.139 


0.139 


Mean velocity 


U 




m s~ 1 


4.31 


8.57 


16.9 


3.12 


5.93 


11.3 


5.59 


11.5 


22.9 


Mean energy dissipation 


(e) = 15v((d !c v) 2 )/2 




m 2 s" 3 


0.141 


0.975 


4.42 


0.244 


2.05 


12.6 


0.379 


2.60 


15.4 


Kolmogorov velocity 


UK = (^<£>) 1/4 




m s~ 1 


0.0375 


0.0610 


0.0891 


0.0428 


0.0736 


0.116 


0.0479 


0.0776 


0.121 


rms u fluctuation 


<M 2 > 1/2 




m s~ 1 


0.236 


0.475 


0.862 


0.552 


1.18 


2.37 


0.743 


1.56 


3.08 


rms v fluctuation 


(v 2 ) 1/2 




m s~ 1 


0.231 


0.463 


0.837 


0.464 


0.975 


1.96 


0.661 


1.36 


2.71 


Skcwncss of u 


<« 3 )/(« 2 } 3 / 2 






+0.07 


+0.05 


+0.06 


-0.22 


-0.12 


-0.10 


-0.03 


-0.04 


-0.04 


Skcwness of v 


{v 3 )/(v 2 ) 3 ' 2 






-0.01 


-0.01 


-0.01 


+0.01 


+0.01 


-0.01 


-0.01 


+0.01 


+0.01 


Kurtosis of u 


<m 4 )/(m 2 } 2 -3 






-0.00 


+0.03 


+0.02 


-0.31 


-0.34 


-0.31 


-0.40 


-0.40 


-0.42 


Kurtosis of v 


(v 4 )/(v 2 ) 2 - 3 






— 0.01 


—0.00 


+0.00 


+0.06 


+0.03 


+0.05 


+0.03 


+0.05 


+0.06 


Kolmogorov length 


r, = (v 3 /(e)) 1/4 




cm 


0.0376 


0.0233 


0.0160 


0.0322 


0.0194 


0.0123 


0.0290 


0.0179 


0.0115 


Taylor microscalc 


A = [2(v 2 )/((d x v) 2 )] 


1/2 


cm 


0.896 


0.684 


0.583 


1.35 


0.996 


0.806 


1.55 


1.21 


0.997 


Correlation length of u 


i« = /" (u(x + r)u 


x))dr/{u 2 ) 


cm 


16.4 


17.0 


18.2 


49.0 


42.4 


43.0 


130. 


128. 


128. 


Correlation length of v 


L v = f£°(v(x + r)v( 


x))dr/(v 2 ) 


cm 


4.13 


4.07 


4.55 


6.94 


6.14 


5.68 


10.3 


10.2 


10.5 


Correlation length of u 2 


L u2 [sec Eq. (1)] 




cm 


4.79 


4.91 


5.40 


18.0 


15.0 


15.0 


21.0 


21.7 


21.9 


Correlation length of v 2 


L„ 2 [sec Eq. (1)] 




cm 


2.42 


2.49 


2.65 


8.33 


7.22 


7.34 


12.8 


12.9 


14.3 


Reynolds number 


Ro A = X(v 2 ) 1/2 /u 






147 


223 


341 


454 


679 


1103 


738 


1183 


1944 


U-u 2 correlation at 4L w 2 


[sec Eq. (9)] 






+0.03 


+0.01 


+0.03 


-0.27 


-0.18 


-0.16 


-0.05 


-0.07 


-0.07 


U-v 2 correlation at 4L^2 


[sec Eq. (9)] 






-0.00 


-0.01 


-0.01 


-0.08 


-0.04 


-0.02 


-0.04 


-0.04 


-0.04 



The maximum is at E = ((N-1)T. In the limit N -+ oo, 
the distribution P(E) becomes Gaussian in accordance 
with the central limit theorem. 3 ' 4 ' 6, 17 

To determine the value of Q, we assume universality of 
P{E) at N 1. Such universality originates in steps of 
the energy transfer. They could be related to interactions 
between the adjacent subsegments. Let us first consider a 
special case where the N subsegments are not correlated 
at all but exhibit the same energy distribution that is for 
the square of a Gaussian random variable. The resultant 
P(E) at any N is the % 2 distribution with N degrees of 
freedom, 17 which corresponds to Eq. (7) with 

C=i (8) 

Then, also at TV ^> 1 in other general cases, the univer- 
sality ensures the same value for (. Since ( = 1/2 yields 
(E) = NT/2 [Eq. (5c)], it might be possible to reformu- 
late the large-scale fluctuations in the classical statistical 
mechanics where (E) = NT/2 holds as the law of energy 
equipartition among N degrees of freedom. 4-6 

IV. CONFIRMATION BY EXPERIMENTS 

The theoretical distribution of E [Eq. (7)] leads to the 
distribution of v^, which is confirmed with experimental 
data of grid turbulence (Gl, G2, and G3), boundary layer 
(Bl, B2, andB3), and jet (Jl, J2, and J3). While Bl, B3, 
and J2 were used in our past work, 9 the others are used 



here for the first time. Their conditions and turbulence 
parameters are summarized in Table I. 

The experiments were conducted under stationary con- 
ditions in a wind tunnel, which had a test section of the 
size of 18 x 3 x 2 m 3 . At a position where the turbulence 
was fully developed, we measured the streamwise veloc- 
ity U + u(t wt ) and the spanwise velocity v(i w t)- Here U 
is the average while u(t wt ) and v(t w t) are temporal fluc- 
tuations. They were converted into spatial fluctuations 
of the longitudinal velocity u(x) and of the lateral veloc- 
ity v(x) by using Taylor's hypothesis, x — —Ut W f For 
further details of the experiments, sec Appendix. 
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FIG. 2. Two-point correlation {[v 2 (x + r) ~ (v 2 )][v 2 (x) - (v 2 )]) 
normalized by its value at r = as a function of r/rj in grid 
turbulence Gl (circles), boundary layer B2 (triangles), and 
jet J3 (squares). The arrows indicate r — 4L V 2. 
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FIG. 3. (Color) Probability distribution of v R /{v R ) at N = 
R/4L V 2 = 10, 30, and 100 in grid turbulence Gl, G2, and 
G3 (black symbols), boundary layer Bl, B2, and B3 (green 
symbols), and jet Jl, J2, and J3 (red symbols). The solid and 
the dotted lines are the theoretical predictions for £ = 1/2 and 
1. 



Figure 2 shows examples of the two-point correlation 
of the local energy of the lateral velocity v 2 , which is used 
to obtain the subsegment length 4L V 2 [Eq. (1)]. The cor- 
relation appears to be almost negligible above the scale of 
4L„2 (arrows) as assumed in our formalism. For a range 
of N — i?/4L„2, we calculate the coarse-grained energy 
v R in each segment with length R [Eq. (2)]. That of the 
longitudinal velocity u is to be studied in Sec. VI. 

Having the length AL V 2 ~ 0.1-0.6 m (Table I), the sub- 
segments of v(x) are local enough to represent local re- 
gions of stationary turbulence that did exist in the wind 
tunnel. The adjacent regions were interacting. We have 
connected these subsegments to make up segments with 
any length R. Although the turbulence in the wind tun- 
nel was limited up to the length of its test section, no 
inconsistency arises because each of the segments is used 
not as a single motion with length R but as a collection 
of N distinct motions with length AL V 2 . In fact, the sta- 
tistical mechanics allows us to make up a canonical sys- 
tem by collecting subsystems that might be even isolated 
from one another, 6 i.e., our subsegments, if they are at 
the same temperature, i.e., T = (v 2 )/y/( [Eq. (5a)]. The 
segment is also homogeneous in the sense that its subseg- 



ments all obey the same statistical law. Since the total 
number of the subsegments is as large as 10 5 -10 6 in the 
individual data records (Table I), the resulting statistics 
are expected to be reliable. 

The mean streamwise velocity Ur in each segment with 
length R is not constant at the value of the mean veloc- 
ity for the entire data, U = (Ur). 18 To confirm that such 
fluctuations of Ur do not affect our study of v 2 R , we cal- 
culate their correlation 

((Ur-(Ur))(v r -(v r ))) (q] 

((Ur- ( Ur )W/ 2 ((vI- K» 2 ) 1/2 ' 1 ' 

which is surely negligible at the scale of the subsegment 
length, R = 4L V 2 (Table I). At around this scale, ((Ur — 
(Ur)) 2 ) is 70% of (u 2 ). 

Figure 3 shows the probability distribution of v R /(v R ) 
at N = R/4L V 2 = 10, 30, and 100. The solid and the dot- 
ted lines are the theoretical predictions of Eq. (7) via Eq. 
(5a) for ( — 1/2 and 1, which depend on N alone. With 
an increase in N, the distribution becomes narrower, but 
it remains wide enough to imply the significance of the 
fluctuations. 7-10 The experiments agree with one another 
and with the theory for ( = 1/2 [Eq. (8)]. 




N 

FIG. 4. (Color) Moments ((v 2 R - (v 2 R )) n ) normalized by {v 2 R ) n 
for n = 2, 3, and 4 as a function of N = R/4L V 2. We also 
show the skewness - (w|>) 3 )/((«| - <"!}) 2 ) 3/2 and the 
kurtosis {(v 2 R - (v R )) 4 ) / ({v 2 R - (v 2 R )) 2 ) 2 - 3, as long as the 
value is reliable. The symbols are the same as in Fig. 3. In 
the top panel, the solid line matches the dotted line. 
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Figure 4 shows ((v R — (v R )) n ) normalized by (v R ) n for 
n = 2, 3, and 4 as a function of N — R/AL V 2. Also shown 
are the skewness ((vl- (v 2 R )) 3 ) / ((v 2 R - (v 2 R )) 2 ) 3/2 and the 
kurtosis {{v 2 R -{v 2 R )) l )/{{v 2 R -(v 2 R )) 2 ) 2 -'i. Theoretically, 
Eq. (7) yields 

((E-(E)) 2 ) =(NT 2 , (10a) 
((E- (E)) 3 ) = 2(NT 3 , (10b) 
((E - (E)) 4 } = (3( 2 N 2 + 6(N)T 4 , (10c) 

with (E) = (NT [Eq. (5c)]. They are used with Eq. (5a) 
to obtain the relations of v R for £ = 1/2 and 1 (solid 
and dotted lines), which depend on TV alone. Again at 
N > 10 1 , the experiments agree with one another and 
with the theory for C = 1/2 [Eq. (8)]. 

Thus, we have reproduced the experiments of the sta- 
tionary homogeneous turbulence. Negligible, if any, are 
undesired variations in conditions of the wind tunnel and 
of the measurement devices. They do not explain the ob- 
served agreement among the experiments, for which v 2 R 
and R have been normalized by (v 2 ) and 4L„2, i.e., char- 
acteristics of the turbulence. We have thereby confirmed 
our thermostatistical formalism as well as its assumptions 
that Rv R is additive and v 2 R has a universal distribution 
at R >• L v 2 . 



V. COMPLETE FORMALISM AND ITS IMPLICATION 

To complete the thermostatistical formalism, we deter- 
mine T/T in Eq. (6). This quantity is dimensionless and 
characterizes the turbulence in the subsegments. Judg- 
ing from T = (v 2 )/y/( [Eq. (5a)], it is natural to equate 
T/Tq with the square of the Reynolds number Re^ for 
those subsegments with length 4L„2 : 



Re 



4L 



l6L 2 v2 (v 2 ) _ ^R 2 T 
v 2 N 2 v 2 



(11) 



Here v is the kinematic viscosity. For £ = 1/2 [Eq. (8)], 
the partition function Z in Eq. (6c) becomes 



Z = Ref L 



(. 



R 2 T 



2„2 



\y/2N'V 



N/2 



(12a) 



The Hclmholtz free energy (F) in Eq. (6b) becomes 

R 2 T 



NT 

(F) — — —— In , _ 
The entropy (S) in Eq. (6a) becomes 



(S) 



(dim _NL( R 2 T 



\ dT J 



\V2N 2 v 2 



(12b) 



(12c) 



Being equivalent to \n(eRe 2 L ) N ^ 2 , the entropy (S) is large 
if the Reynolds number Re4L is high. Lastly, the resis- 
tance force (/) is obtained as 



f d(F) \ _ NT 
V dR ) T ~ ~R' 



(12d) 



This is analogous to a force originating in the Reynolds 
stress, d Xj (viVj) for the velocity Vi in a direction Xi. The 
reason is NT / R oc (v 2 )/L v 2, where L v 2 serves as the 
scale for a significant variation of v 2 . 

Our formalism of Eq. (12) agrees with the thermody- 
namics. While T and (/) are intensive, (S), R, and (F) 
are extensive, cx N. Through the Legendre transforma- 
tion, (F) yields other thermodynamic potentials, 5,6 e.g., 
the Gibbs free energy (G) = (F) + (f)R as a function of 
Tand(/): 



(G) 



NT 



2 -In 



2^3 



V2v 2 {f) 2 



(13) 



with (S) = — (dr(G))(fy These potentials are totally dif- 
ferentiate and hence reproduce the Maxwell relations. 5,6 
For example, from [d R (d T (F)) R }T = [d T {d R (F}) T ] R in 
Eq. (12), we have 



d(S) \ 
dR ) T 



[dim 

\dTj R 



(14a) 



To describe an equilibrium state, the potentials also re- 
produce the so-called thermodynamic inequalities. 5,6 For 
example, between C R = T(d T (S)) R from (F) in Eq. (12) 
and C(/) = T(9t(S')){/) from (G) in Eq. (13), we have a 
relation analogous to a well-known inequality 5,6 between 
heat capacities at constant pressure and at constant vol- 
ume: 

3iV N 
C (f) = ^>C R =->0. (14b) 

This agreement with the thermodynamics implies that 
our formalism is surely thermostatistical. The conclusion 
remains the same even if £ ^ 1/2, for which we only have 
to consider Z = Re^^ and so on. 



VI. CONCLUDING DISCUSSION 

For stationary and homogeneous turbulence, we have 
studied large-scale fluctuations of the coarse-grained en- 
ergy of the lateral velocity v 2 R [Eq. (2)]. They have been 
described in a thermostatistical formalism, which has the 
same mathematical structure as used for the statistical 
mechanics of equilibrium systems with many degrees of 
freedom. By using an analogy between the fluctuations of 
v 2 R [Eq. (3b)] and the thermal fluctuations of the energy 
E [Eq. (4)], we have obtained a correspondence between 
v R and E [Eq. (5)]. The resultant formalism reproduces 
the distribution of v R observed at N = R/4L V 2 > 10 1 
in Figs. 3 and 4 [Eqs. (7) and (8)]. Therefore, through 
the large-scale fluctuations, turbulence is related to the 
statistical mechanics. 

The thermostatistical formalism of Onsager for a class 
of two-dimensional turbulence is well known. 3,13,19 We 
have demonstrated that such a formalism also exists at 
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large scales of the usual three-dimensional turbulence, al- 
though we have used a canonical ensemble while Onsager 
used a microcanonical ensemble 4-6 at constant E. 

To construct the formalism, we have assumed that Rv 2 R 
is additive at R ^> L v 2 [Eqs. (3b) and (5a)], by assuming 
that any n-point spatial correlation of v 2 is negligible at 
r ^> L v 2 . We have also assumed that the distribution of 
v R is universal at R 3> L v 2 [Eq. (8)], by assuming that the 
dependence on the flow configuration is lost after many 
random steps of the energy transfer to the large scales. 
These assumptions have been confirmed along with the 
formalism. Especially in Figs. 3 and 4, we have observed 
the universal distribution of v R . 

However, the additivity and the universality might not 
be exact. The spatial correlations of v 2 might not be ex- 
actly negligible in flows such as those known to have or- 
ganized motions far above the scale of the energy supply 
L, 20,21 where the additivity might be lost in some manner 
specific to the flow configuration. Since the effective de- 
grees of freedom of the segment might be less than TV, the 
effective value of C might be less than 1/2 if we recall the 
discussion leading to Eq. (8). Thus, our formalism might 
not be exact. Nevertheless, it is at least a good approx- 
imation because we have observed no disagreement with 
the experiments. Similar discussions exist about possible 
effects of the flow configuration on otherwise universal 
motions at small scales. 3,10 ~ 12 ' 22 

The additivity and the universality have not been con- 
firmed for the coarse-grained energy of the longitudinal 
velocity u 2 R . Figure 5 shows that our formalism does not 
reproduce the experiments of some of the flows. Their 
skewness and kurtosis, among others, are greater than 
those for £ = 1/2 (solid lines), implying that their effec- 
tive values of £ are significantly less than 1/2. These are 
significant versions of the above mentioned feature. Pre- 
sumably, along the one-dimensional cut of an incompress- 
ible fluid, longitudinal distortions propagate to larger dis- 
tances than lateral distortions. The spatial correlations 
of the local energy u 2 do not always decay fast enough 
to become negligible at the scales studied here. Still at 
the larger scales, there remains a possibility to reproduce 
experiments of all the flows. 

Our formalism does not hold at R < 476^2 , but an ap- 
proximation is available. For example, if the distribution 
of v is Gaussian, we have ((v R — (v R )) 2 ) —> 2(v 2 ) 2 in the 
limit R —¥ 0. Then, Eq. (3a) is approximated at any R 
as {{v R - (v R )) 2 ) = 4:L v 2(v 2 ) 2 /(R + 2L v 2). By comparing 
this with Eq. (4), we obtain correspondences as in Eq. (5) 
and a formalism as in Eqs. (6)-(8) and (12). 

There is also an application of our formulation to large- 
scale fluctuations other than those of v R , if they are in 
an equilibrium as well as have additivity and are thereby 
analogous to thermal fluctuations in the statistical me- 
chanics. We only have to write the mean square in the 
form of Eq. (3b) and compare it with Eq. (4) . The formu- 
lation could be constrained by some additional feature, 
e.g., universality in Eq. (8). This is the case not only 
for one-dimensional data as studied here but also for any 




N 

FIG. 5. (Color) Same as in Fig. 4 but for u 2 R as a function 
of N = R/AL U 2 . For the kurtosis, we only show the values of 
Gl, G2, G3, and Jl. These are the most reliable experiments 
where the correlation of Ur with u R at R — 4L U 2 is negligible 
(Table I). Between the values of G1-G3 and Jl, scattered are 
the values of the other experiments. 

data of the higher dimension. Examples are expected to 
be found in a variety of fluctuations, far beyond those of 
turbulence. 
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Appendix: DETAILS OF EXPERIMENTS 

The experiments were conducted in a wind tunnel of 
the Meteorological Research Institute. We adopt coordi- 
nates x w t, Uwt, and z w t in the streamwise, spanwise, and 
floor-normal directions. The origin x w t = y w t — z wt — 
m is on the floor center at the upstream end of the 
test section of the wind tunnel. Its size was 5x w t = 18 m, 
Sy w t — 3 m, and 5z w t = 2 m. The cross section 5y w t x Sz wt 
was the same upstream to x w t — —4 m. 

The wind tunnel had an air conditioner. If needed, we 
used this conditioner to constrain the variation of the air 
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temperature. The resultant variation was ±1 °C at most 
in each experiment, where the kinematic viscosity v is 
assumed to have been constant. 

To simultaneously measure U + u and v, we used a hot- 
wire anemometer. The anemometer was composed of a 
constant temperature system and a crossed-wire probe. 
The wires were made of platinum-plated tungsten, 5 fim 
in diameter, 1.25 mm in sensing length, 1mm in separa- 
tion, oriented at ±45° to the streamwise direction, and 
280 °C in temperature. 

For the grid turbulence, we placed a grid at x w t = — 2 
m across the flow passage to the test section of the wind 
tunnel. The grid had two layers of uniformly spaced rods, 
with axes in the two layers at right angles. The cross 
section of the rod was 0.04 x 0.04m 2 . The spacing of the 
rod axes was 0.20 m. We set the incoming flow velocity to 
be 4 (Gl), 8 (G2), or 16ms- 1 (G3). The measurement 
position was on the tunnel axis, y w t = 0m and z wt 
1 .01) mi (Table I). 

For the boundary layer, roughness blocks were placed 
over the entire floor of the test section. The block size was 
Sx w t — 0.06 m, Sy w t = 0.21m, and <5z w t = 0.11m. The 
spacing of the block centers was <5a; w t = Sy wt — 0.50 m. 
We set the incoming flow velocity to be 4 (Bl), 8 (B2), or 
16ms _1 (B3). The measurement position was in the log- 
law sublayer at x wt = +12.5 m and y wt = 0m (Table I), 
where the boundary layer had the displacement thickness 
of 0.2 m and the 99% velocity thickness of 0.8 m. 

For the jet, we placed a contraction nozzle. Its exit 
was at x wt = — 2 m and was rectangular with the size of 
Sy wt — 2.1m and <5z wt = 1.4m. The center was on the 
tunnel axis. We set the flow velocity at the nozzle exit to 
be 8 (Jl), 16 (J2), or 33ms- 1 (J3). The measurement 
position was at x wt = +15.5 m, y wt = m, and z wt — 
0.40 m. 

These measurement positions were determined so that 
the skewness (v 3 ) /(i> 2 ) 3 / 2 and the kurtosis (v 4 ) / (v 2 ) 2 — 3 
were close to the Gaussian value of (Table I). It en- 
sures that the turbulence was fully developed and various 
eddies filled the space randomly and independently 23,24 
Not always close to the Gaussian value were (w 3 )/(u 2 ) 3 / 2 
and (u 4 )/(u 2 ) 2 — 3 (Table I). They are sensitive to spe- 
cific features of the energy-containing eddies that depend 
on the grid, the roughness, or the nozzle. 

The signal of the anemometer was linearized, low-pass 
filtered, and then digitally sampled. We set the sampling 
frequency as high as possible (Table I), on the condition 
that high-frequency noise was not significant in the power 
spectrum. The filter cutoff was at one-half of the sam- 
pling frequency. We obtained a long record of 1 x 10 8 or 
4 x 10 8 data in each of the experiments (Table I). 

The sampled signal is proportional to the flow velocity, 
through the calibration coefficient that depends on the 
condition of the anemometer and thereby varied slowly 
in time. For individual segments of each data record, the 



length of which is fixed for the record and ranges from 
4 x 10 6 to 2 x 10 7 data, we determined the values of the 
coefficient so as to have the same U value. The coefficient 
within any of these segments is estimated to have varied 
by ±1% at most. 
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